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Abstract 

We calculate the universal flavor-singlet radiative QED corrections to unpolarized lep- 
ton scattering applicable to general differential scattering cross sections, involving 
charged fermions or photons in initial or final states. The radiators are derived to 
0((a ln(Q 2 /m^)) 5 ) in analytic form. Numerical illustrations are given. 



1 Introduction 



QED corrections to integral and differential cross sections of light charged lepton-anti-lepton 
scattering or deeply inelastic lepton-nucleon scattering turn out to be quite large in some kine- 
matic regions [1-4]. This applies in particular to the Bremsstrahlung contributions due to signifi- 
cant shifts in the kinematics of the underlying differential scattering cross sections. The universal 
corrections can be grouped into flavor non-singlet and flavor singlet contributions. In the orders 
0((e 2 aL) k ), with a the fine structure constant, L = \n(Q 2 /m 2 ), Q 2 the typical virtuality of the 
process and to/, e/ the fermion mass and charge, respectively, the non-singlet contributions stem 
from the leading order anomalous dimension in QED, Pff(x,Q 2 ). The diagonal elements of the 
singlet anomalous dimension matrix contain a 5(1 — ^-distribution and are distribution- valued 
due to (...) + -distributions. Therefore the numerical Mellin-inversion does badly converge in the 
region of x J$ 1. Analytic representations are required to high enough order k in the fine struc- 
ture constant a to cover all universal effects for the energy ranges probed at present day colliders 
and those to be built in the foreseeable future. This applies to high luminosity experiments at 
future linear colliders [5] and as well to the search for rare reactions at LHC. The second order 
universal corrections for various processes are known for a long time [6-9]. The 3rd order cor- 
rections were given in [10] 1 for the flavor non-singlet and in [12] also for the singlet case. Later 
the 5th order non-singlet corrections were given in [13] and recalculated in [14] and [15,16] 2 . In 
Ref. [15] a very compact form was given for the non-singlet contributions, which are the same 
for polarized and unpolarized scattering. There also the polarized singlet contributions were 
calculated to 0((aL) 5 ). In the present paper the unpolarized singlet evolution kernels are calcu- 
lated to 0((aL) 5 ) which supplements earlier investigations for the flavor non-singlet kernel [15]. 
A second class of universal QED corrections was treated previously in the non-singlet [17] and 
polarized singlet case [15]. It concerns the leading order small x resummations of 0((a\n 2 (x)) k ). 
These resummations are based on corresponding resummations in QCD [18]. In the unpolarized 
(pure) singlet case the leading order small-x QCD resummations [19] result from the non-abelian 
gluon coupling, which is absent in QED. Therefore a transformation of the respective kernels is 
not possible in this case. 

The paper is organized as follows. The general framework to derive the radiators in 0((aL) k ) 
is outlined in section 2. In section 3 the leading order singlet radiators are calculated to 0((aL) 5 ). 
In section 4 numerical illustrations are given and section 5 contains the conclusions. An appendix 
lists useful convolution relations which were needed in the calculation and are of use in other 
QED and QCD calculations. 

2 The Solution of Singlet Evolution Equations 

The universal QED corrections 0((aL) k ) for general values of the collinear radiation momentum 
fraction x can be expressed solving the singlet evolution equations starting at a low scale Qq. This 
scale may be identified with a typical charged lepton mass toj squared. The running coupling 
constant a(Q 2 ) = a(Q 2 )/(An) obeys the evolution equation 




(1) 



1 For an application to the Z-pe&k see [11]. 

2 The results of [13] and [15] agree but partly disagree in the 5th order with [14]. 
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where (3k denote the expansion coefficients of the QED /^-function in the MS-scheme, @ = 
-A/3N f ,p 1 = -4:N f ,p 2 = 2N f + (44/9) JVj etc. [20] for N f active charged lepton species. At 
leading order the solution 

a(m 2 f ) 

with L = ln(Q 2 / mj) is obtained. 

The universal radiators are found as solutions of the leading order QED renormalization 
group equations associated to the collinear singularities. The following QED-distributions are 
of relevance : 

D^ s (a, x) = D^(a,x) — D^(a,x) (3) 
= D f (a,x) + D f (a,x) (4) 
D 11 (a,x) = D 22 (a,x) (5) 
D lf (a ) x) = D 21 (a,x) (6) 
Df~f(a,x) = D 12 (a,x) (7) 

The non-singlet distribution D^ s (a, x) was dealt with in [13-15] before. Here the Dij(a, x) denote 
the respective matrix element of the singlet radiator Ds(a,x) given below. 
The singlet radiator functions at the scale Qq are 

D(Ql)(x)= D(a„) = l 6(1 -x) (8) 

since both the charged leptons and the photon are considered to be asymptotically stable parti- 
cles. The evolution equations read 



da a J2T=o ah Pk 

1 



(3 a 

The Mellin convolution is defined by 



P (x) + a ( P x (x) - |- P (x)j + 0(a 2 ) 



D s (a,x). (9) 



A(x)<8iB(x)= I dxi I dx 2 5(x - x 1 x 2 )A(x 1 )B(x 2 ) . (10) 
Jo Jo 

The Mellin transform 

M\f(x)](N)= f dxx N - 1 f{x) (11) 



Jo 

diagonalizes the convolution (10) to 

M[A(x) ® B(x)](N) = M[A(x)](N) ■ M[B(x)](N) . (12) 

In some of the radiators (...)+ distributions emerge, which are defined relative to the set of 
smooth test functions <j)(x) with compact support by 



I 



dx[F(x)] + <j>(x) = [ dxF(x) [<f>(x) - 0(1)] . (13) 
o ' Jo 



Eq. (9) may be solved easiest in Mellin space as a matrix-valued ordinary differential equation 
to all orders, see e.g. [21,22]. The leading order solution reads 

D s ,o{a,x) = [exp(- R (x)ln(a/a ))®\ <g> D s (a Q ,x) = E {a,a ,x) <8> D s {a , x) , (14) 

where R = Pq/Po and E denotes the leading order singlet evolution operator. We use the 
short-hand notation 

[/(^))®] = E^r^®, fc =i»(x), (is) 

fc=0 

with <S>f =1 the /c-fold convolution. The singlet solution (14) to fcth order in a(Q 2 )L therefore 
requires to calculate fc-fold convolutions of the leading order matrix of splitting functions. 

The method described above can be extended to sub-leading logarithmic contributions, i.e. 
terms of 0(a 2 L) etc. These contributions contain pro cess- dependent parts being described by 
Wilson coefficients in inclusive situations or the corresponding semi-inclusive quantities. These 
corrections are neither universal nor independent of the measurement chosen for the kinematic 
variables. Examples are the 0(a 2 L) corrections for the initial state radiation in e + e~ annihilation 
[6], the 0(a 2 L) initial and final state radiation corrections to deeply inelastic scattering [23] and 
the the 0(a 2 ln(m M /m e )) corrections to the electron spectrum in muon decay [24]. 



3 The Leading Order Solution to 0[(cxL) 5 ] 

In the following we derive the solution for the unpolarized singlet QED evolution kernels up to 
0((aL) 5 ). The singlet evolution equation is solved in the running coupling a(Q 2 ). However, 
one may re-parameterize the representation and express the evolution kernel directly in terms 
of do = a{m 2 ) by 

E (a,x) = 1 5(l-x)+ P (x)a L+ Q P«( x ) + P (^ (a L) 2 (16) 
IpW(x) + |p«(x) + ^Po(x)) (a L) 3 
+ (1 P ( J\x) + i P^(x) + I P«(x) + P P (s)) (a Ly 

+ (m p ° ){x) + 1 pi ° ){x) + 1 p?)(x) + !? p ° ){x) + ii Po{x) ) (a ° L)5 • 



Here the matrices P^ are 



Pfti)=®<"P (i) , (17) 



i.e. ®V>P (x) = [P ® P ] (x) etc. The corresponding expressions for the (1, l)-components of 
P^\x) are given relative to the non-singlet components P^g(x), see [15]. The singlet-matrices 
are 

P ° {X) - [ P^\x) P^\x) J ' (18) 
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with components P^ (x) given below in Eqs. (19-39). They were calculated using the convolution 

formulae of Appendix A and relations given in [15,25] before. The projections P-p describe the 
splitting of a fermion into a fermion (1,1), of a photon into a fermion, positron, respectively (1,2), 
a fermion into a photon (2,1), and a photon into a photon (2,2) in kth order in the renormalized 
coupling constant. 

The leading order QED splitting functions can be obtained identifying Tr — Cp — 1 and 
Ca = in the QCD splitting functions [26,27] in accordance with the gauge group U(l). 

1st order terms: 



p£\ x ) = P^)+P^(x) = ^^--l-x + l5(l-x) (19) 

PllUx) = (20) 

P$(x) = 2[x 2 + (l-x) 2 ] = Ax 2 -Ax + 2 (21) 

P£\x) = 1 + ^~ X)2 =x-2+ 2 - (22) 

X X 

Pg\x) = —5(1- x) (23) 

The 5(1 — x) distribution in (23) emerges due to momentum conservation for the photon mo- 
mentum Jq dxx [Pi2(x) + ^22(2^)] = 1- 

The 2nd order terms are : 



pS\ x ) 



3(2). 
12 1 



x) 



P\msi x ) + -furs 
ln(l — x 



,(2) 



X 



8 



1 — x 
6 



-4(1 + a;) ln(l - x) + \n(x) 



7 + 7x- 



1-x 



1-x) 



8 9 8 
3 - Sx - -x 2 + — + 



!j-4C(2) 



5(1 -x) 



7 20 



4(1 - 2x + 2x 2 ) ln(l -x)- 2(1 - 2x + Ax 2 ) \n(x) - - + — x - -x 2 



x-2 + 



x 



P£\x) = 21n(l-x) 

p£\x) 

They were derived in [7] and various other places before. 



3 3 

\ 1 z \ 10 7 4 
+ (2-x) \n(x) + — - -x 
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3a; 



4(1 + x) \n(x) + 2-2x--x 2 + — + -5(1 - x) 

3 3x 9 



(24) 
(25) 
(26) 
(27) 



3rd order terms: 



Pii( x ) — Pims( x ) + Piips( x ) 



,(3) 



(3) 



24 



ln 2 (l 



x 



1-x 
-ln(l -x) 



- 12(1 + x) ln 2 (l -x) 
32 32 



36 



22 + 14a: + ^x 2 - — 

3 6X 



-241n(x) 



ln(l -x) 

1-x 
ln(l - x) 
l-x 



5 



x 



+34(1 + x) ln(x) ln(l - x) + ln 2 (x) 



ln(x) 



83 59 



32 



18 



1 x H x — 

6 6 3 1-x 



4 15 15 

i - x ~ Y~Y X 



+ 22(l + x)Li 2 (l-x) 



+ 



'27 

y - 24C(2) 
'27 



(1-x). 



- - 18C(2) + 16C(3) 



317 155 16 2 16 in/i 
-- + -* + Y * 2 -- + 12(l + x K (2) 

8(1 -x) 



8(1 - 2x + 2x 2 ) ln 2 (l - x) - -(5 - 16x + Ax 2 ) ln(l - x) 

o 

-8(1 - 2a; + 4x 2 ) ln(x) ln(l - x) - (3 - 6x - 8x 2 ) ln 2 (x) 



ln(x) 



32 26 _ 2 

1 £ — lbx 

3 3 



16x 2 Li 2 (l - x) 



521 491 232 2 32 . rt rt 

+ x x 2 + 8(1 - 2x + 2x 2 )((2) 

18 9 9 9x V AV ; 



x-2 + 



x 



ln 2 (l -x) - 



5x - 16 + 



x 



ln(l -x) 



3 1 
+4(2 - x) ln(x) ln(l - x) + -(2 - x) ln 2 (x) + -(26a: - 19 

2 3 



16 



) ln(x) 



+ -Li 2 (l -x) -4 

x 



P 2 ( 2 3) (x) = 41n(l-x) 



x-2 + - 

x 

4 „ 41 
l-x--x 2 + -- 
3 3 x 



(28) 



(29) 



491 521 16 2 116 

C(2) + x + —x 2 

sv ; 18 36 9 9x 

+ 8(1 + x) ln(x) ln(l - x) - 2(1 + x) ln 2 (x) 



(30) 



31 



32 



--(4 + x - 4x 2 ) ln(x) + 8(1 + x)Li 2 (l - x) (1 — x) H 



2 1 
x 

X 



-—S(l-x) 
27 K J 



(31) 



4th order terms: 



Here and for the 5th order terms we refer to the expressions P^s* ( x ) gi ven i n R- e f • [15] for brevity. 



P^x) = P£i(x)+ 48(1 + x)ln(x)ln 2 (l -x) - 24(1 + x) ln 2 (x) ln(l - x) 



,(4) 



— (1 + x) ln 3 (x) + 8 
3 



3 - 3x - 4x 2 + - 

x 



ln 2 (l -x) 



—8 ln(x) ln(l — x 
-4 



"4 


14 o 2 _ 




25 






—x - 8x 2 


+ ln 2 (x) 




+ -x — 16x 2 


3 ~ 




~y 





73 _ 73 

¥ ~ Y x 



16 



-x 



16 

9 ~ ' 9x 
569 1445 128 



ln(l - x) 

64 



-y + — x + T" x2 + fx + 48C(2) + m2)x 



+96(1 + x) ln(l - x)Li 2 (l - x) + 8 



55 ,2 
- + -x + 4x 2 + - 

3 3 x 



ln(x) 
4 



Li 2 (l-x) 



745 745 

-96(1 + x)Li 3 (l - x) + 48(1 + x)S lj2 (l - x) - — + — x + 

lo lo 



224 2 224 
~ X ~ ~9x 
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P$(x) = 16(1 + x) ln(x) ln 2 (l - x) - 8(1 + x) ln 2 (x) ln(l - x) - 2(1 + x) ln 3 (x) 

16 

ln 2 (l - x) - — (2 - x - 4x 2 ) ln(x) ln(l - x) 



3 - 3x - Ax 2 + 



x 



C(2) 



pg) = 16(l-2x + 2x 2 )ln 3 (l-x) - 24(1 - 2x + 4x 2 ) ln(x) ln 2 (l - x) 



-2 ln 2 (x) ln(l - x)(5 - lOx - 24x 2 ) + ln 3 (x) 
"13 44 



7 14 



16 



x x 

3 3 3 



-41n 2 (l -x) 



+ ln 2 (x) 



80 



4 ln(x) ln(l — x) 



32 35 



-x — 16x 



79 14 

1 X X" 

6 3 3 



ln(l -x) 



959 1978 



9 



+ ln(x) 



1505 614 



-x + 



832 



320 



-x 



128 
~9x~ 



+ 48C(2) - 96C(2)x + 96C(2)x 2 



18 



9 



x + —x 2 + 24C(2) - 48C(2)x + 96C(2)x' 



-96x 2 ln(l - x)Li 2 (l - x) - 44(1 - 2x) ln(x)Li 2 (l - x) 



-4Li 2 (l - x) 



40 2 

15 — 3x — — x 

O 



+ 96x 2 Li 3 (l - x) 



-4(13 - 26x + 16x 2 )Si i2 (l -x) + 



18065 5947 



+4 



13 44 8 2 



108 

C(2) + 32(l-2x + 2x 2 )C(3) 



560 2 

x H x 

27 9 



128 

27x 



P 2 ? = 8 



x-2 + - 

x 



ln 3 (l - x) + 12(2 - x) ln(x) ln 2 (l - x) 



+5(2 - x) ln 2 (x) ln(l - x) - 7 -(2 - x) ln 3 (x) + \ ln 2 (l - x) 



44 - 13a: - - 



x 



-2 ln(x) ln(l — x 

~989 959 



53 58 32 1' 

3 3 3 x 



+ ln(l -x) 



64 



9 



-x H x — 

18 9 3x 



+ ln 2 (x) 
160 



16 101 



12 



-x 



48 



+ 48C(2)-24C(2)x--C(2) 



x 



682 413 64 2 64 

+ x x 2 + 24C 2 + 12C 2 x 

9 36 9 9x SV ; SV ; 



+ ln(x) 

+ — ln(l - x)Li 2 (l - x) + 22(2 - x) ln(x)Li 2 (l - x) 

x 



-2Li 2 (l - x) 



3 - 15x + 



40 

3x 



48 

--Li 3 (1-x +2 

x 



26 - 13x + 



x 



(32) 



(33) 



Si )2 (l-a;) 



5947 18065 64 „ 280 

I x x H 

54 216 27 9x 



44 - 13x 



x 



C(2) 



-16 



x 



X 



C(3) 



(34) 



+8 



1 - x - -x 2 + 



7 



--(19 - 23a; + 16a; 2 ) ln 2 (x) - 41n(l - x) 



16 , 
9 - 9a; - — x 2 + 
9 



16 

9x 



— ln(x) 



?H + ^ + ^ + !i + 16C(2) + 16C(2)l 



8. 



+32(1 + x) ln(l - x)Li 2 (l -x)- -Li 2 (l - a;) 
-32(1 + x)Li 3 (l - x) + 16(1 + ar)Si >2 (l - a;) 



1 + x - 4x 2 

x 



1331 1331 



18 



608 2 608 
+ ^ X + ^7 X ~27x 



l - x ~r + Tx 



C(2) 



+—5(1 - a;) 
81 V ; 



(35) 



5th order terms: 



P[f(x) = Pg£(x) + 128(1 + x) ln(x) ln 3 (l - x) - 96(1 + x) ln 2 (x) ln 2 (l - x) 



,(5) 



+-(l + x) ln 4 (x) + 64 
3 



1 — X X 

3 



4 

3a; 



ln 3 (l - a;) 



-32(1 - 5a; - 8a; 2 ) ln(x) ln 2 (l - x) - 8(7 + x + 16a; 2 ) ln 2 (x) ln(l - x) 



+ - [10 - 5a; + 16a; 2 ] \n 3 (x) - 64 
9 



x 2 1 
6 — 6a; 1 

3 3a; 



ln 2 (l -x) 



-32 \n{x) \n(l-x 
+ 161n 2 (a;) 
-4 ln(l -x) 



97 313 



9 



9 



-x 2 + — + 12C(2) + 12C(2)x 



CO 1 oo 

^ + ^ + 2a ; 2 + 6C(2) + 6C(2)a; 



331 331 



9 



7T X ~ 

9 



1024 



-x z + 



1024 
27a 7 



64 



+48C(2) - 48C(2)a; - 64C(2)a; 2 + -((2) 



+4 ln(x) 



1990 997 

+ —r^X 



320 



64 



9 '°° + 27a; 



27 27 

+8C(2) - 40C(2)a; - 64C(2)a; 2 + 64((3) + 64C(3)a;] 
+384(1 + x) ln 2 (l - x)Li 2 (l - x) - 96(1 + x) ln 2 (a;)Li2(l - x) 



+ 128 



l+x + 2x z + - 

x 



ln(l — a;)Li 2 (l — x) 



-144(1 - x) ln(a;)Li 2 (l - x) - 768(1 + x) ln(l - x)Li 3 (l - x) 
+384(1 + x) ln(l - ar)Si, 2 (l - x) - 256(1 + x) ln(a;)S li2 (l - x) 



-32 



-128 



205 205 n 2 2 
— + —x + 2x 2 + - + 12C 2 + 12C(2 )x 
9 9 x 



1 + x + 2x 2 + 



x 



Li 3 (l -x) - 16 



Li 2 (l-x) 

7-15x-^x 2 -^ 
3 3x 



Si, 2 (l-x) 



+768(1 + x)Li 4 (l - x) - 192(1 + x)Si, 3 (l - x) - 384(1 + x)S 2)2 (l - x) 



8 



+ 27 
+ 128 



10127- 10127x- 



1 — x x H 



2456 



-x z + 



2456 
3x 



+ 64 



6 — 6x — —x + 



3x 



C(2) 



3a; 



C(3) 



(36) 



P$(x) = 32(1 -2x + 2x 2 )ln 4 (l -x) - 64(1 - 2x + 4x 2 ) ln(x) ln 3 (l - x) 

g 

-24(1 - 2x - 8x 2 ) ln 2 (x) ln 2 (l - x) + -(5 - lOx - 16a; 2 ) ln 3 (x) ln(l - x) 



+ ln 4 (x) 



12-r+r 



64 

^-(2-7x + x 2 )ln 3 (l-x) 



-16(8 + llx - 12a; 2 ) ln(x) ln 2 (l - x) + -(26 + 23a; - 60a; 2 ) ln 2 (x) ln(l - x) 

3 



+ ln 3 (x) 
-4 ln 2 (l -x 



, x 224 2 

5 1 x 2 

3 9 



649 1478 728 2 32 
"9 9 _;r + ~ ;r "3^ 



+48C(2) - 96C(2)a; + 96C(2)a; 2 ] 



+41n(a;) ln(l - x) 



1027 866 



9 



x + 

9 



1328 



9 



x 1 + 48C(2) - 96C(2)a; + 192((2)a; 2 



+ ln 2 (x 
+41n(l -x) 



451 1664 



1072 



-x — 



9 9 9 

"6064 7721 2056 
-x + 



x 2 + 24((2) - 48C(2)a; - 192((2)a; 2 



+ \n(x) 



27 27 
8861 24511 



x 2 - 



128 



+ 32C(2)-112C(2)a;+16C(2)a; 2 



54 



+ 



27 



x - 416a; 2 



27 27a; 

+64C(3) - 128C(3)a; + 128C(3)x 2 ] 

256 



27a; 



+ 128C(2) + 176C(2)a; - 192C(2)a; 2 
■128C(3) + 256C(3)a; - 512C(3)a; 2 ] 



-384x 2 ln 2 (l - ar)Li 2 (l - x) + 288(2a; - 1) ln(x) ln(l - x)Li 2 (l - x) 

+ 16(1 - 2a; + 4a; 2 ) ln 2 (a;)Li 2 (l - x) - 32(12 - 3a; - 10a; 2 ) ln(l - x)Li 2 (l - x) 
32 

+ — (1 - 5a; - 12a; 2 ) ln(a;)Li 2 (l - x) + 768a; 2 ln(l - x)Li 3 (l - x) 

-288(2a; - 1) ln(a;)Li3(l - x) - 32(11 - 22a; + 16a; 2 ) ln(l - x)Si, 2 (l - x) 
-16(5 - lOx - 16a; 2 ) ln(a;)Si i2 (l - x) 



+8 



100 ,16 . . o 

21 + 34a; + —x 2 + — + 48( (2 )x 2 
3 3a; 



U 2 (l-x) 



+32(12 - 3a; - 10x 2 )Li 3 (l - x) - 16 



10 + 5a; x 



Si )2 (l-a;) 



-768a; 2 Li4(l - x) + 32(1 - 2x + 2x 2 )Li 2 (l - x) 

-32(7 - 14x - 6x 2 )S lj3 (l -x)+ 32(7 - 14x + 8x 2 )S 2>2 (l - x) 



+ 



322519 240997 



648 
128 



324 



-x + 



23776 



81 



-x 



1088 4 
127a 7 + 3 



649 1478 



-x + 



728 



-x 



32 

x 



(2 - 7x + x 2 )C(3) + 48(1 - 2x + 2x 2 )C(4) 



C(2) 
(37) 



9 



-16 



2 -x - 



x 



ln 4 (l - x) + 32(2 - x) ln(x) ln 3 (l - x) 



20 5 

+ 12(2 - x) ln 2 (x) ln 2 (l - x) (2 - x) ln 3 (x) ln(l - x) - — (2 - x) ln 4 (x) 

3 24 

32 T 1 1 8 

+ — 7-2x-- ln 3 (l-x)-8 17 - 16a: + - 

3 a; x 



ln(x) ln 2 (l -x) 



+^(41 - 46a;) ln 2 (x) ln(l - x) + -(x - 47) ln 3 (x) 

O 



+ ln 2 (l - a;) 
-ln(x) ln(l -x) 
-In 2 (a; 



2956 1298 



64 



1456 



9 



-x H x — 

9 3 9a; 



192 



+ 192C(2) - 96C(2)x C(2) 

Ob 
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The radiators can be expressed through these functions and polynomials thereof as well as 
rational functions in x. 

The universal radiator functions Eq. (14) can now be attached to the respective initial- or 
final-state radiating light charged fermion or photon lines of any differential scattering cross sec- 
tion to account for the respective leading order QED corrections. These radiators generalize the 
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radiators due to soft-photon exponentiation, valid for Dns((i(Q 2 ),x) [29] in the region x — > 1, 
to general values of x and all collinear transitions possible. The numerical effect of the respec- 
tive radiator depends on the change of the subsystem kinematics of the differential scattering 
cross section, which usually differs for initial and final state radiation and due to the type of 
leg encountered. This kinematics has to be worked out for the respective process accordingly. 
Moreover, the radiative corrections may strongly depend on the way the kinematic variables of 
the process are measured. In case of deeply inelastic scattering investigations of these aspects 
were performed in [3,8,9,28]. For similar considerations for e + e~ annihilation see e.g. Ref. [2]. 

The radiative correction due to the radiator D aia2 (a(Q 2 ), x) for a differential cross section 
reads 



(41) 



b\=b\...bi=b\ 



Here the / kinematic variables which determine the differential cross section are b r \ r=1 . Their 



rescaled value under changing the momentum p, 



final state radiation, b r 
corresponding Jacobian 



-> z • p ai resp. p ai — > p ai /z for initial or 
, is bounded by Zq 1 for hard radiation. J ai (b r \ l r=1 , z) denotes the 



J ai (brl =1 ,z) 



dh/dh 



dk/db, 



dk/dk 



(42) 



and d l a a2 /dbi...db[ is the subsystem differential cross section for which the line of type a± is being 
replaced by a line of type a 2 . Eq. (41) may be generalized to the case of more universal radiators 
correspondingly, requiring additional rescaling of variables. 

In the above we assumed, that the radiator functions describe collinear radiation along outer 
fermions or photons. However, in various applications also internal, nearly collinear situations 
may occur. One example is the so-called 3rd [30] or Compton peak [2, 3, 31] in deeply inelastic 
scattering. Here a photon being originally virtual contributes near to its mass shell in the radia- 
tive correction, which gives rise to a factorizing collinear process. The universal contributions to 
these processes can be obtained from radiators as well. 

Finally we would like to add a remark on small x resummations. In QCD the leading order 
corrections stem from 
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and the function 71, (iV, a s ) obeys 

l L (N,a s ) = 



C A a s 



7T(N 
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21+1 



(45) 



i=i 



with a s = a s /(4ir) the strong coupling constant and (k Riemann's (^-function. The transition to 
QED, Ca — ► 0,CV — > 1 trivializes (43,44) except for the lowest order term in a(Q 2 ) in P^°(x) 
which is already known from the fixed order terms above. Yet all anomalous dimensions do 
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receive 1/x terms in higher orders, which contribute, as well-known [32], in the abelian limit 
starting from next-to-leading order : 
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However, these terms do not stem from the resummation [19] but are of different origin. Their 
pole strength is of 0(a 2 / (N — 1)), which is larger than for the poles resulting from 0(a 2 ln 2 (a;)) 
terms. For unpolarized QED radiators systematic resummations of the leading small x terms 
were not carried out yet. At leading order in a(Q 2 ) only P 7 /(x) oc 1/x, while at next-to- 
leading order (46-49) all terms contain a singular contribution. It would be worthwhile to derive 
resummations of these terms in the future. 

4 Numerical Results 

The singlet contributions to the universal radiator functions, summing the leading logarithmic 
corrections up to 0((aL) 5 ), are shown in Figure 1, as a function of the momentum fraction x for 
different values of Q = \[Q 2 in the case of mj = m e . The corrections in case of other charged 
fermions have to be rescaled accordingly in L = ln(Q 2 /m 2 ). In case of fractionally charged 
leptons, a has to be replaced by e 2 ao. Here Dff denotes the pure singlet part of the fermion 
radiator to which the non-singlet contribution has to be added, cf. [15]. The diagonal radiator 
vanishes at leading order 0(aL), while the radiator D 22 contributes at x — 1 only due to 
momentum conservation. For x < 1 the (pure) singlet (PS')-diagonal radiators contribute with 
0((aL) 2 ) only. In this order the PS'-diagonal terms are identical in the region x < 1 and differ 
by the 5(1 — ^-distribution in D 22 . At higher orders both diagonal elements receive different 
corrections. This explains the relative smallness of the radiators D\f and D 22 compared to D\ 2 
and D 2 i. Yet the diagonal radiators grow oc 1/x as x — > 0. This growth is even more pronounced 
in D 2 i, which contains a 1/x term already at 0(aL), while D12 receives those terms at 0((aL) 3 ) 
only and therefore shows a moderately varying profile in x. The QED scaling violations shown 
in Figure 1 are of moderate size, comparing scales from Q = 10 GeV to Q — 1 TeV which is due 
to the smallness of the fine structure constant and its weak running. At x — 0.1 the radiators 
-Dfi and -D22 reach about 1% and grow further towards smaller values of x. D 12 takes values 
between 3 and 8%. _D 2 i is largest and reaches 50% at x — 0.1. The radiators D\{ 3 , D 2 \ and D 22 
vanish towards x — > 1, while D 12 approaches finite values. 

Figure 2 compares the size of the first order contribution to D u , D± 2 and D 2 \ with the respec- 
tive contributions up to 0((aL) 5 ), which can be regarded numerically as the total contribution 
for the values of Q chosen. For Du the first order contribution is much smaller than the total 
contribution in the region of small values x due to the steep rise of the pure singlet component. 
The higher order contributions to D\ 2 are small at medium values of x, and amount to a —5% 
correction in the small x region and a +10% correction at large values of x for Q = 10 GeV. 
The higher order contribution to D 2 \ range between 1% to 10%. 
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Figure 3 shows the impact of the 5th order term w.r.t. the first four orders. The effect of 
the 5th order term amounts to O(10~ 5 ) for D u , D 12 and D 2 \, while the effect in the case of D 22 
is one order of magnitude larger at large value of x. In either case, we confirm that the singlet 
radiator to the 5th order has a very high accuracy, which is sufficient to represent the universal 
part of QED corrections, relevant in high precision measurements in both high-energy charged 
lepton anti-lepton collision, charged lepton-nucleon collisions and photon collisions at future 
electron-positron linear colliders with a possible Giga-Z option, electron-photon and photon- 
photon colliders and future muon colliders. These are the reactions to which the corresponding 
radiators are contributing. The reactions do also contribute to the precision measurements of 
QCD scaling violations [3, 28, 33] in deeply inelastic scattering as universal QED corrections. 
Likewise they are important for rare initial states at high energy hadron colliders such as LHC 
and would contribute in lepton and photon initiated processes there, such as single leptoquark 
production [34]. 

The radiator functions calculated above are made available in form of a FORTRAN-program 
which can be obtained form the authors upon request. 

5 Conclusions 

The collinear logarithms in QED can be resummed due to the renormalization group equations 
for mass factorization. Unlike the case in QCD the collinear logarithms are finite due to the 
possibility to define the coupling constant asymptotically, i.e. in the limit of vanishing scales. 
The associated logarithms are well defined since photons and leptons are non-confined and may 
be regarded as stable or long-lived states. The leading order corrections 0((aL) k ) are universal. 
The respective radiators resum the radiative corrections which only depend on the type of particle 
transition % — > j. As shown, sufficient numerical stability of O(10~ 4 " -5 ) is reached evaluating the 
radiators to 0((aL) 5 ) for scales as large as Q < 1 TeV. The radiators are presented in rr-space 
and can be applied directly to the respective multiply differential scattering cross sections to 
describe the universal contribution due to initial and final state radiation off the different outer 
legs contributing to the respective scattering process involving charged fermions and photons. 
In the small x region the leading order radiators receive contributions cx 1/x, with an onset in 
different orders in aL, which leads to larger corrections in this kinematic region. A systematic 
resummation of the particular small x contributions, unlike the case for the non-singlet and 
polarized singlet corrections, is not known yet. The radiators derived can easily be adopted for 
experimental analysis and simulation programs. 
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6 Figures 
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Figure 1: The singlet radiators as a function of x and Q in %. D^f denotes the pure singlet 
part of D u . 
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Figure 2: Relative contribution of the first order singlet radiators D\- in all terms to 0((aL) 5 ). 
Here Dn denotes the sum of the non-singlet contributions -D NS , with soft exponentiation beyond 
0((aL) 5 ), and the pure singlet contribution -Dff. 
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Figure 3: Relative contribution of the singlet radiators -D^- up to the 4th order in aL if compared 
to all terms to the 5th order as a function of x and Q, with D u = D^ +PS . 
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7 Appendix A: Mellin Convolutions 



In this appendix we list the convolutions of functions up to weight 5 required in the present 
calculation in addition to those given in Ref. [15,25]. Some of the integrals require to use 
Mellin transforms and algebraic relations between the finite harmonic sums [35,36]. They were 
calculated recursively in explicit form and may be of general interest for other higher order 
calculations in QED and QCD. 
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